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1.  Evaluation  of  Analytic  Electronic  Wavefunctions  and  Properties  for 
Atoms  and  Molecules:  Theoretical  Development 

The  theoretical  formulation,  as  well  as  computer  programs  for  the 
evaluation  of  Hartree-Fock  wavefunctions  and  various  properties  of  atomic  and 
molecular  systems,  has  been  completed.  Extensive  calculations,  carried  to  the  limit 
of  the  Hartree-Fock  independent  particle  model  demonstrated  the  need  for  the  explicit 
inclusion  of  electron-electron  correlation  into  the  electronic  wavefunction,  if  various 
atomic  and  molecular  properties  were  to  be  computed  reliably.  To  this  end  the  multi¬ 
configuration  self-consistent-field  (MC-SCF)  model  was  developed.  In  this  model  the 
electronic  wavefunction  is  represented  by  o  linear  combination  of  several  Slater 
determinants,  rather  than  only  one,  as  in  the  Hartree-Fock  model,  and  oil  the  orbitals, 
as  well  os  the  configuration  expansion  coefficients,  are  determined  simultaneously 
using  the  variational  priixiple.  The  development  of  this  model  proceeded  in  two 
stages:  fint  the  method  was  restricted  to  pair  replacement,  i.e.  the  shell  occupation 
numben  in  the  different  Slater  determinants  of  the  wavefunction  had  to  differ  by  even 
numbers.  Later  this  restriction  was  removed  and  a  completely  general  theory  was 
developed.  Gsnsiderable  effort  was  expended  towards  the  development  of  stable 
numerical  procedures  for  the  solution  of  the  strongly  coupled  MC-SCF  equations. 

This  development  was  completed  successfully.  Computer  programs  for  the  calculation 
of  MC-SCF  wavefunctions  for  atoms  and  linear  molecules  were  constructed  and  tested 
extensively.  The  theory  for  the  calculation  of  excited  states  of  some  symmetry  as 
lower  lying  states  was  developed.  This  has  been  incorporated  into  the  MC-SCF 
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computer  programs,  thus  the  calculation  of  excited  states,  even  if  there  are  lower 
states  of  same  symmetry,  can  be  performed  with  equal  ease  and  accuracy  as  the 
calculation  of  ground  states. 

Methods  for  the  optimization  of  the  nonlinear  parameters,  the  orbital 
exponents,  in  analytical  wavefunctions  have  been  improved,  thus  allowing  a  more 
rapid  accurate  calculation  of  atomic  and  molecular  wavefunctions.  Also  the  calcu¬ 
lation  of  wavefunctions  for  linear  molecules  has  been  made  significantly  more 
efficient  through  a  completely  new  analysis  of  the  required  molecular  integrals. 

This  analysis  was  carried  out  in  collaboration  with  Drs.  McLean  and  Yosimine  of 
IBM  Research,  and  has  lead  to  a  molecular  integrals  program,  which  Is  significantly 
faster  and  more  accurate  than  similar  programs  avjilable  heretofore. 

2.  Calculation  of  Electronic  Wavefunctions  and  Properties 
of  Atoms  and  Molecules 

Using  the  available  Hartree-Fock  computer  programs,  the  calculation 
of  analytical  wavefunctions  close  to  the  Hartree-Fock  limit  has  been  completed  for 
fint  and  second  row  diatomic  hydrides.  For  these  molecules  potential  curves  and 
various  properties  were  computed.  The  characteristic  of  the  chemical  binding  and 
the  electronic  charge  distribution  has  been  analyzed  extensively.  Electron  affinities 
for  CH,  NH,  OH  and  SiH,  PH  and  SH  have  been  computed  in  the  Hartree-Fock 
approximation,  and  corrected  empirically  for  the  correlation  error  inherent  in  the 
Hartree-Fock  approximation.  In  the  case  of  polyatomic  systems,  wavefunctions  for 
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CH^  and  OH^  were  obtained  using  the  one  center  expansion  method  and  the 
Hartree-Fock  approximation. 

A  careful  study  of  the  2  and  4  electron  atoms  and  ions  has  been 
undertaken,  in  order  to  assess  the  value  of  the  MC-SCF  model  and  its  efficiency  in 
accounting  for  electron-electron  correlation.  Preliminary  calculations,  using  the 
MC-SCF  model,  have  been  performed  for  several  low  lying  states  of  LiH  and  BeO. 

These  calculations  are  still  preliminary.  However,  indications  are  that  in  the  case 
of  LiH  potential  curves  for  several  low  lying  states  can  be  obtained  with  an  accuracy 
of  a  few  100  cm  ^  without  excessive  computational  effort. 

3.  Calculation  of  Transition  Moments 

Computer  programs  for  the  calculation  of  absolute  line  strength  of 
electronic  and  vibration-rotational  specho  of  diatomic  molecules  have  been  developed. 
These  programs  require  as  input  information  the  molecular  electronic  energy  and  wave 
function  at  various  internuclear  distances.  With  this  the  rotation-vibration  wave 
functions  and  line  intensities,  f-numbers,  are  calculated. 

To  date  the  vibrational  transition  probabilities  of  the  electronic 
ground  states  of  BeO,  MgO,  CaO  and  SrO  have  been  calculated  using  HF  electronic 
wave  functions  and  potential  curves  as  well  as  experimental  RKR  potential  curves.  The 
detailed  results  of  these  calculations  have  been  written  up  for  publication. 
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In  addition,  the  cnicuiation  transition  moments,  in  dipole  length 
and  momentum  representation  of 
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is  in  progress.  A  report  of  the  results  obtained,  together  with  a  critical  review  of  the 
reliability  of  measured  and  computed  transition  probabilities,  has  been  completed. 

4.  Numerical  SCF  for  Atoms 

An  extensive  and  thorough  study  was  carried  out  re-examining  the 
classical  approach  of  Hartree  and  Fock  to  calculate  atomic  wove  functions.  In  the 
course  of  this  work,  many  refinements  to  the  theory  were  discovered  which  Influence 
numerical  computations  considerably,  especially  for  heavier  atoms. 

The  generally  accepted  formulas  for  the  orbitals  In  the  regions  r  -•  o 
and  r  »  were  found  to  be  In  error;  the  correct  formulas  were  derived,  and  they 
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yielded  some  surprising  results.  For  example,  it  was  found  that  for  r  -*  <»  all 
orbitals  ultimately  decay  with  a  decay  constant  appropriate  to  the  weakest  bound 
electron. 

New  formulas  and  techniques  were  derived  which  facilitated  the 
adjustment  of  orbital  energies,  and  also  for  matching  the  inner  and  outer  portions 
of  the  wave  functions. 

Additional  work  is  necessary  before  a  computer  program  can  be 
completed,  based  on  this  analysis.  This  is  expected  to  be  accomplished  in  the  near 
future.  Calculations  of  much  greater  accuracy  and  reliability  than  carried  out  to 
date  will  then  be  possible. 

5.  Variational  Treatment  of  Electron-atom  Scattering  and 
Atomic  Resonance  States 

Two  complimentary  efforts  were  undertaken  in  this  area.  One  involved 
the  numerical  determination  of  orbital  and  scattering  wave  functions,  the  other  on 
expansion  technique  similar  to  the  LCAO-SCF  procedure.  Both  efforts  ore  based  on 
o  wave  function  of  great  flexibility,  in  order  to  overcome  the  inherent  limitations 
in  close  coupling  ond  related  methods  commonly  used.  Much  attention  was  paid  to 
avoid  the  numerical  difficulties  encountered  by  singular  behavior  of  the  mathematical 
formalism,  which  have  so  often  been  encountered  in  calculations  of  this  type.  Papers 
describing  these  efforts  will  be  forthcoming  in  the  near  future. 


7 


6.  Relativistic  Gilculation  on  Atoms 

The  commonly  used  approaches  to  calculate  relativistic  wave  functions 

2 

for  atoms  were  carefully  examined.  It  was  found  that  up  to  order  (v/c)  the 
relativistic  corrections  are  best  calculated  by  the  straightforward  generalization  of 
the  Pauli  approximation  to  the  many-elactron  case.  Hence  all  relativistic  corrections— 
to  this  order  of  accuracy— can  be  determined  as  expectation  values  using  non-relativistic 
Schrodinger  wave  functions  and  appropriate  operators.  It  Is  also  shown  that  this  formalism 
can  be  applied  to  wave  functions  of  great  generality;  in  particular,  MC-SCF  wave 
function  can  be  used  with  ease  In  order  to  calculate  relativistic  corrections  with  great 
accuracy. 

Calculations  of  level  splittings  for  several  atoms  with  low  Z  have  been 
carried  out,  and  are  continuing.  A  paper  describing  the  theory  and  the  corresponding 
calculations  is  nearing  completion. 

7.  Atomic  and  Molecular  Polarizabilities 

Pblarizabillties  of  atoms  had  been  calculated  successfully  In  this 
loboratory  using  the  coupled  Hartree-Fock  formalism.  Present  Investlgatlore  concen¬ 
trate  on  using  wave  functions  of  much  greater  generality.  Very  occurate  polarizabilities 
have  been  calculated  for  a  few  light  atoms,  and  vrark  on  other  atoms  and  the  N2 
molecule  is  In  progress. 
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In  connecHon  with  this  work,  several  Interesting  theoretical  results 
have  been  derived,  using  Hellman-Fe/nman — ,  hypervirlal  — ,  and  other  well-known 
theorems. 

8.  Accurate  Wave  Functions  for  the  H2  Molecule 

Further  work  on  the  H2  molecule.  In  particular  for  excited  states  of 
Interest,  was  carried  out  by  Drs.  Kolos  and  Wolniewicz.  AH  this  work  has  now  been 
published,  and  need  not  be  summarized  seporately. 
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